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Abstract 

We demonstrate a novel method of introducing point defects (mono and 
di-vacancies) in a confined mono-layer colloidal crystal by manipulating in- 
dividual particles with optical tweezers. Digital video microscopy is used to 
study defect dynamics in real space and time. We analyze the topological 
arrangements of the particles in the defect core and establish their connection 
with the energetics of the system. It is found that thermal fluctuations excite 
point defects into dislocation multipole configurations. We extract the dislo- 
cation pair potential at near field, where cores overlap and linear elasticity is 
not applicable. 

PACS number(s): 82.70.Dd, 05.40.Jc, 61.72.Ji, 61.72.Bb, 61.72.Ff 
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Colloidal crystals [|T|, ordered self- assembled structures of (sub)micron spheres, provide 
a model system for the study of basic problems in condensed matter physics. In particular, 
confined colloidal crystals have been used to explore statistical physics in two dimensions 
[||,|^. Although colloids in 2D have been studied for the last 20 years, several areas still 
remain relatively unexplored. The simplest structural defects in a crystal, namely vacancies 
and interstitials, fall in this category. Point defects are of considerable interest, since they 
play a dominant role in real materials' properties. Also, they are predicted to proliferate 
in more exotic systems, such as the ground state of a Wigner crystal or the supersolid 
phase [0 of the Abrikosov lattice in type II superconductors. 

Colloidal crystals also offer new approaches in synthesizing materials with novel proper- 
ties and applications like photonic-band-gap materials , optical switches , and chemical 
sensors p. Understanding the behavior of defects in colloidal crystals, and inventing tech- 
niques for manipulating their dynamics can have an immediate impact on the fabrication of 
nanostructured materials. 

In this work we demonstrate the use of optical tweezers 0, a powerful tool that has 
promoted research in a wide range of fields, to artificially introduce isolated point defects 
in otherwise structurally perfect two-dimensional colloidal crystals. For the first time, such 
defects are studied by video microscopy in real space and time, providing new insight into 
their microscopic dynamics. In this letter the energetics of the defects are discussed. Results 
on the diflFusion of the defects are presented in separate publication [|10[ . 

The colloidal crystals were prepared with a ^1% Vol. aqueous suspension of 0.3 /xm 
diameter negatively charged polystyrene-sulfate micro spheres (Duke Scientific, polydisper- 
sity The ionic strength of the suspension was reduced by ion exchange to a few fiM. 

Under these conditions the screening length is of order a few hundred nm and the particles 
crystallize due to strong electrostatic interaction. Confining the suspension between two 
fused silica substrates separated by ^2 /xm suppresses the vertical motion of the spheres by 
negative charge that develops at the silica- water interface and creates a single layer colloidal 
crystal with a lattice constant a = 1.1/xm [pT[ . 



Trapping a particle with optical tweezers and dragging it from its lattice site creates 
isolated point defects in an otherwise perfect two-dimensional crystal The optical 



tweezers were formed by focusing a beam of an Ar+ laser (Coherent INNOVA 90, A =514nm) 
through the same objective (Zeiss Plan Neofiuar, lOOx, oil-immersion, NA=1.3) used for 
imaging. The position of the beam was controlled by a set of mirrors mounted on an XY 
galvanometer scanner set (Cambridge Technology). Point defects are created by dragging a 
particle away from its lattice point at about 50 /xm/sec, making it an interstitial and leaving 
a vacant site behind [|TB|. To overcome the restoring force from the crystal and to be able 



to drag the particle faster than the relaxation time of the lattice, 100 mW of laser power 
50 mW going through the objective lens) was used. 

The dynamics of such single, isolated point defects were studied in real time under 100 x 
magnification and recorded on videotape (Sony SVO-9500MD recorder) using a monochrome 
CCD video-camera (Sony SSC-M370). The field of view covers (50 x 40)a^ Typically, we 
processed 60 frames/sec and tracked a single defect for 20-40 seconds, which is the time it 
takes to diflFuse away from the field of view. Each individual frame was acquired on a PC and 



processed using a particle-tracking algorithm [14|. The raw data consists of the positions of 



the particles in every frame which are linked into trajectories in time. 
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Real-space imaging facilitates detailed study of the structure of the defect core. The fol- 
lowing questions naturally arise: what are the possible configurations in which a defect can 
exist? How can these be identified? What are their symmetries? Are there any geometrical- 
topological constraints that must be satisfied? In the case of strong interactions, na < 5.9, 
where is the screening length, numerical studies of mono-vacancies and intersti- 

tials have revealed that the system spontaneously deforms into configurations with lower 
symmetry than that of the lattice. This prediction is verified for our system, where hia ^ 3 
(Fig. 0A-D). For di- vacancies, various stable configurations with well defined symmetry can 
also be identified (Fig. gE-H). 

The interactions in the system are short ranged (screened Coulomb, V{r) oc exp{—nr)/r)^ 
therefore the nearest neighbor bonds have geometrical as well as physical significance, since 
they identify the pairs of mutually interacting particles. The creation of a point defect 
involves breaking of bonds and creation of mis-coordinated particles in the core region. 
The configuration of the defects can be characterized by the arrangement of these broken 
bonds and mis-coordinated particles. We analyzed the topological structure of the defects 
by performing a geometrical triangulation of the positions of the particles, finding the pairs 
of nearest neighbors, and measuring the coordination number rii for every particle. Our 
observations indicate that the following constraints are satisfied: (i)the mis-coordinated 
particles are never observed isolated but rather appear as distinguishable pairs, triplets, 
etc. of nearest neighbors, (ii)the mean coordination number of every pair, triplet, etc. is 
always equal to six, (iii)we can define vectors fiij starting from a particle i with < 6 and 
terminating to a nearest neighbor j with rij > 6. The number of vectors originating from 
i is 6 — and the number of vectors terminating on j is — 6. For every configuration. 

The finite temperature causes the system to vibrate around every local energy mini- 
mum C, exploring the volume of phase space Qc in which the energy is 0{kBT) above the 
minimum energy Ec^ and occasionally getting enough energy to jump to a nearby local mini- 
mum. The contribution to the entropy of every configuration is roughly Sc = klog{Qc)^ and 
the diflFerent configurations occur with relative probabilities Pc oc exp{—Fc/kBT)^ where 
Fc = Ec — TSc is the free energy of a configuration. As long as the system remains around 
a local energy minimum, the distortions of the lattice are elastic and the topological ar- 
rangement of the particles does not change. Therefore, every energy minimum has not only 
a characteristic symmetry but also a certain topological configuration of the defect core. 
Using the topology rather than the symmetry as a criterion for identifying the different 
configurations C, we measure their relative occurence probabilities Pc- The free energy 
diflFerences, estimated from the measured Pc's, are found to be of order 0.1-1 fc^T. 

Point defects are "topologically neutral" since they have zero Burgers vector. However, 
one could think of creating a point defect by inserting extra rows of particles that do not 
terminate on each other. The extra inserted rows are edge dislocations whose Burgers 
vectors add up to zero. Therefore, a point defect contains embryonic pairs, triplets etc. of 
dislocations, in the same way that a dislocation contains an embryonic pair of a seven and 



a five- fold disclination [|1^]. Since the system is at finite temperature, one would expect the 
defect to fluctuate between these possible configurations of dislocations. 

The signature of an isolated edge dislocation in a two-dimensional hexagonal lattice is 
a core with a pair of neighbor 5-fold and 7- fold coordinated particles. The extra row of 
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particles terminates on the 5-fold particle. The Burgers vector 6 = a is in 1-1 correspon- 
dence with the vector fiij defined earlier in the text, b connects the 5-fold coordinated 
particle with its first nearest neighbor, moving counter-clockwise after the 7-fold neighbor 
(see Fig.|^). Since b ^ Hij and Z)<ij> ^ij ^ 0, one can view the various configurations of the 
defect as X] 6 = dislocation multipoles. Tracking the time evolution of the system, we see 
that the vast majority of the observed configurations correspond to a dislocation pair or a 
triplet, with infrequent appearance of higher order multipoles. Most of the time the dislo- 
cations comprising the defect are closely bound in the configurations previously identified 
(Fig.|). Occasionally however configurations are observed in which the dislocations appear 
to dissociate and recombine (Fig.|D. This fits well with a picture in which one sees the dislo- 
cations as point particles interacting with an attractive potential. The bound configurations 
are minima of the energy and the system gets random kicks from the finite temperature 
background getting exited into states in which the dislocations can be separated. 

Focusing on the case in which the point defects appear as a dislocation pair and assuming 
an interaction V{r) between the two dislocations, we would expect Boltzmann statistics for 
the probability of observing the two dislocations a certain distance apart, namely: P{r) oc 
exp{—V{r)/kBT). In addition, particle conservation dictates that dislocations can only glide 
parallel to their Burgers vector, so r = Tc+r^, with Vc ^ 3a(2a) for di (mono) vacancy and Vg \\ 
b. We measured P(r), r = |f], for the mono and di- vacancies, using a few thousand snapshots 
of the system. The separation r of the two dislocations is identified as the separation between 
the two five-fold coordinated particles where the extra rows of particles terminate. Our 
results (Fig.^ show a rather rapid decrease in P(r) as r goes beyond a couple of lattice 
constants, together with a modulation, more clearly seen in the case of the di-vacancy. 
Since V{r) oc —logP{r)^ the interaction of the two dislocations turns out to increase with 
r, with an average slope 2.9 ± 0.4 and 2.0 ± 0.2kBT/a for the case of mono and di-vacancy 



respectively. The modulation corresponds to a Peierls-like energy barrier [[l^] that comes 
from the discreteness of the lattice. A final point is that the core region of dislocations has 
a typical size of a few a. Therefore, our measurements of V{r) are in a regime where the 
cores of the two dislocations overlap and linear elasticity theory is not applicable. 

In summary, we demonstrate a novel way of introducing point defects in colloidal crystals 
through manipulation of individual particles with optical tweezers. Using digital video 
microscopy, we identified the topological configurations of the particles around the defect and 
showed how they are connected with the energetic configurations found by numerical work 
in [|1^ . Topological analysis of the possible configurations reveals that the point defects can 
appear as dislocation multipoles with zero total Burgers vector. The individual dislocations 
are bound by an attractive interaction. Di- vacancies in particular are observed to dissociate 
into two dislocations that can separate several lattice spacings before recombination. 

We acknowledge helpful discussions with Prof. S.C. Ying. This work was supported by 
NSF (grant DMR-9804083), the Petroleum Research Fund, and the Research Corporation. 
XSL acknowledges the support of the A. P. Sloan Foundation. 
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FIGURES 



FIG. 1. (LeftjExperimental setup (not to scale): The sample cell consists of a 1/2" diameter 
quartz disk(Q) and a quartz coverslip(CS) glued together. The distance between the two surfaces 
was controlled by a patterned thin 2fim) polymer film (S) (Dow Chemical Co., CYCLOTENE), 
which served both as a spacer and as an adhesive. The cell was connected to an external circulation 
circuit containing a few mL of suspension, driven by a peristaltic pump (P) (VWR) and including 
a conductivity meter (M) (VWR, model 1054, flow-through cell) and a column with mixed-bed 
ion-exchange resin (lEX) (BIO-RAD, AG 501-X8(D)). (C) Indicates a 0.5mm circulation channel 
in contact with the two-dimensional region and (T) a particle trapped with the optical tweezers. 
(T^i^M^Micrograph of an isolated di- vacancy in a two-dimensional colloidal crystal. The 4x4 
diamond in the enlarged central region contains 4 particles in the perfect lattice but only 2 when 
it encloses the core of the di- vacancy. 

FIG. 2. Configurations of mono-vacancy: (A)split (SV)^ (B) symmetric (V3), (C,D) crushed 
(^2a5^26)- Nomenclature adopted from Ref.l5. Top left insets in every figure show configura- 
tion of the vectors fiij defined in the text. Top right insets show the arrangement of miscoordi- 
nated particles in the core of the defect. The relative free energies of the different configurations 
were estimated to be (in kBT): Fsv — 0.86, Fya = 1.10, Fy2a + Fy2h = 1.41. Configurations of 
di-vacancy: (E,F)split (S'De,^)? (C) crushed (1^2)? (D) symmetric (£^3). Free energies in UbT are: 
Fd2 = 0.7, FsD = 1.1, Fd3 = 2.1. 

FIG. 3. Point defects as dislocation multipoles: (A) well-defined di-vacancy configuration. A 
Burgers circuit (solid line) fails to close if it crosses the core of the defect. A Burgers circuit that 
surrounds the defect core without crossing it (dashed line) closes as expected. The arrows indicate 
the Burgers vectors of the "embryonic" dislocation dipole. (B) configuration of the same di-vacancy 
about 2 seconds later. The di-vacancy is not well defined any more and the system resembles a pair 
of dislocations with opposite Burgers vectors, whose gliding lines (dotted lines) are separated by 
^ 3a. When those two dislocations come close to each other, the original di-vacancy is recovered. 

FIG. 4. Pair interaction between dislocations, di(top) and mono-vacancies(bottom). P(r) was 
determined from a few thousand snap-shots of the system by measuring the distribution of the 
separation of the two 5-fold coordinated particles found. V{r)/kBT was estimated from —logP{r). 
The solid curves are spline fits and the straight lines are linear fits to the particular regions of data. 
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